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The wetting properties of solid substrates with customary (i.e., macroscopic) random roughness 
are considered as a function of the microscopic contact angle of the wetting liquid and its partial 
pressure in the surrounding gas phase. Analytic expressions are derived which allow for any given 
lateral correlation function and height distribution of the roughness to calculate the wetting phase 
diagram, the adsorption isotherms, and to locate the percolation transition in the adsorbed liquid 
film. Most features turn out to depend only on a few key parameters of the roughness, which can 
be clearly identified. It is shown that a first order transition in the adsorbed film thickness, which 
we term 'Wenzel prewetting', occurs generically on typical roughness topographies, but is absent on 
purely Gaussian roughness. It is thereby shown that even subtle deviations from Gaussian roughness 
characteristics may be essential for correctly predicting even qualitative aspects of wetting. 
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I. INTRODUCTION 

While the physics of wetting and spreading on ideally 
smooth solid surfaces has meanwhile reached a status of 
mature textbook knowledge, a whole range of wetting 
phenomena on randomly rough substrates are still elu- 
sive. This is particularly annoying, as almost all surfaces 
of practical interest bear considerable roughness, be it 
due to weathering, wear, or on purpose as, e.g., in the 
case of sand-blasted surfaces. Clearly, this has substan- 
tial impact in many situations. For example, a drop of 
liquid deposited on a rough substrate will spread or not, 
depending on the morphology of the liquid film which de- 
velops in the troughs of the roughness. As it accommo- 
dates its free surface to the substrate topography, it may 
percolate across the sample. The drop will then gradu- 
ally spread over the entire sample. On the contrary, if 
the film rather tends to form isolated domains, the drop 
will stay in place. Similarly, the redistribution of liquid 
within a granular pile, such as in humid soil or sand, may 
proceed along the grain surfaces only if the liquid wetting 
film on the grains forms a percolated structure. The mor- 
phology of a liquid water film deposited from humid air 
onto the surface of an electric isolator will strongly affect 
the performance of the latter, for analogous reasons. 

There has been already a lot of research on the adsorp- 
tion of liquids on rough surfaces [1-15], but this was con- 
cerned either with roughness amplitudes as small as the 
(nanometer) range of van der Waals forces [31 El El 12] or 
with rather artificial substrate topographies in the con- 
text of super-hydrophobicity [T0HT31 H"5] . or with com- 
pletely wetting liquids (zero contact angle) [J SI El M- 
The most frequently encountered, customary case, how- 
ever, is characterized by a finite contact angle, and a 
random roughness with typical length scales at least in 
the micron range. In the present paper, we consider sur- 
faces which exhibit a random topography on scales large 
as compared to molecules, and are subject to adsorption 
of a liquid which forms a small but finite contact angle 
with the substrate material. Since we consider macro- 



scopic roughness, we adopt the view that all interfaces 
are infinitely sharp on the length scale of consideration 
(sharp kink approximation [16 ). As the typical length 
scales considered here are still small as compared to the 
capillary length of the liquid (2.7 mm for water), gravity 
will be neglected as regards its effect on the liquid surface 
morphologies to be described. As opposed to earlier stud- 
ies which concentrated on the macroscopic contact angle 
and contact line [H [5] , we will try to derive the wetting 
phase diagram and other characteristics connected to the 
adsorption of a liquid film. 

The first systematic study of wetting on a randomly 
rough substrate at finite contact angle owes to Wenzel 
|17) . He characterized the roughness by a single parame- 
ter, r, which he defined as the ratio of the total substrate 
area divided by the projected area. Obviously, r > 1, 
and r = 1 corresponds to a perfectly smooth surface. 
The free energy which is gained per unit area when the 
rough substrate is covered with a liquid is then given by 
r(7sg — 7si)> where 7,,; and j sg are the solid-liquid and 
solid-gas interfacial tension, respectively. If this is larger 
than the surface tension of the liquid, 7, we expect a 
vanishing macroscopic contact angle, because covering 
the substrate with liquid releases more energy than is re- 
quired for the formation of a free liquid surface of the 
same (projected) area. More specifically, force balance 
at the three-phase contact line yields 

COS Omacro = " — T COS 6 (1) 

7 

for the macroscopic contact angle on the rough surface. 9 
is the microscopic contact angle according to Young and 
Dupre. When the microscopic contact angle is reduced 
to 9w = arccos(l/r), which we will henceforth call Wen- 
zel's angle, 9 macro vanishes, and the substrate is covered 
with an 'infinitely' thick liquid film. As we will see be- 
low, however, there are imprortant ramifications which 
are sensitive to the kind of roughness of the substrate. 
Furthermore, even minute deviations from liquid-vapour 
coexistence, as they are omnipresent in practical situa- 
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tions, unveil a rather complex scenario which goes well 
beyond eq. 



II. PRESENTATION OF THE PROBLEM 

We describe the topography of the rough solid sub- 
strate by f{x,y), where (x, y) is a vector in the plane. 
The (randomly varying) function / is normalized such 
that < / >= 0, where the angular brackets denote aver- 
aging over the entire sample area, iS. It is assumed that 
the substrate is homogeneous and isotropic, in the sense 
that the statistical parameters of / are the same every- 
where on the sample, and independent of rotation of the 
sample about the normal axis of the sample. 

A small amount of liquid deposited on this substrate 
will make an interface with the surrounding gas, which 
can be described by a second function, g(x,y). The sup- 
port of g is the set W C S, which denotes the wetted 
area. Continuity of the liquid surface assures that g = f 
on the boundary of W, i.e., at the three-phase contact 
line, where the solid substrate, the liquid, and the gas 
phase meet. This line will henceforth be denoted by dW- 

Information on g(x,y) can be obtained from the total 
free energy functional of the system, which is given by 



w 



F : / hsi - ls g )V^ + (V/) 2 + 7V 7 ! + (V.9) 2 ] dx dy 

Minimization of F yields two important properties of g. 
First of all, the mean curvature of the liquid surface, 
which can be written as 1181 



Based on the observation that the amplitude of most 
natural roughness is much smaller than its dominant lat- 
eral length scale, we assume for the present study that 



V/ |« 1 



(4) 



which allows for substantial simplifications. Expanding 
then the mean curvature according to eq. ([3]), we obtain 
to first order in Vg 



H w -A<? 
2 y 



(5) 



Similarly, the contact angle with the substrate yields the 
boundary condition 



V(g-f) | wtaafl 



(6) 



on dW, to first order in V/ and Vg. We can now im- 
mediately write down a useful identity concerning these 
quantities. Green's theorem tells us that 



/ nV( 5 -/)d,s= f A( 5 -/)d 2 x 
Jaw Jw 



(7) 



where s is the distance along <9W, n its unit normal vec- 
tor, and x = (x,y). Since g = f on 9W,V(g — /) is 
everywhere perpendicular to dW. Hence eq. ^ may be 
written as nV(g — /) rj tan#, and we can rewrite eq. ^ 
as 



L tan 9 + [ [2H - A/]d 2 x = 
Jw 



(8) 



in which L denotes the length of dW. This equation will 
be the starting point of the discussion to follow. 



2 I v/l + (V.g) 2 



(3) 



assumes the same value everywhere on W. Second, the 
two surfaces described by f(x,y) and g(x,y) make the 
same (Young-Dupre) angle 9 everywhere on dW . This 
reflects the local force balance at the three-phase contact 
line. 

While surface roughness gives rise to substantial con- 
tact angle hysteresis on macroscopic scales, the micro- 
scopic contact angle, 9, is known to be well defined on the 
typical (micrometer to nanometer) scale [19l [20] . Never- 
theless, we should be aware that even on small scales, 
equilibration will take time, be it by transport through 
the gas phase or through an adsorbed layer of molecu- 
lar thickness [91 [TB] (which we disregard in the present 
study). 

The question we shall ask is the following. Given 
the substrate topography, f(x, y), the equilibrium micro- 
scopic contact angle, 9, and the mean curvature of the 
liquid surface, H, what can we predict on the function 
g(x, y) and the shape of the wetted area, W? In particu- 
lar, we shall be interested whether W forms a percolated 
set in the plane. 



III. GAUSSIAN ROUGHNESS 

If we want to exploit eq. QSb , we have to refer to a 
specific roughness function, j(x). Following the over- 
whelming majority of the literature on randomly rough 
surfaces, we will start by considering Gaussian roughness. 
The height distribution is then 



P(f) 



1 



exp{-/ 2 /2C } 



(9) 



where Co is the mean square of /(x). / can be fully 
characterized by its lateral correlation function j^T] [2"T1 
[23H26], 

C(r)=</(x)/(x + r)> (10) 
Below we will make use of its polynomial 



with r =| r 
expansion, 



C(r) = Co 



n=l 



C n r r ' 



(11) 



where Co is the mean square amplitude of the roughness. 
The form of eq. ( 10 1 reflects the isotropy of the roughness 
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A. Distribution functions 

For Gaussian roughness, the joint distributions of / 
with other stochastic quantities can be obtained in a 
straightforward manner from multivariate analysis. As 
it is well known [2TJ [M], the joint distribution of two 
quantities f\ and / 2 is then given by 



and 



exp 



(12) 



where {Mij) is the inverse of the matrix (< fyfj >) and 
Q is the determinant of that matrix [2Tj. For the joint 
probability of / and V/, we find 



K/>v/) 



K/) 
471-a 



exp 



(V/) 2 
4C, 



(13) 



On the side, this directly yields < (V/) 2 >= 4C 2 , from 
which we can conclude that roughness topographies ful- 
filling @ will have C 2 < 0.003. 

For the joint probability of / and A/, we obtain 



p(/,A/) = 



2ttVQ 



exp 



64C 4 / 2 -8C 2 /A/ + C (A/) 2 
2Q 



(14) 

with Q =| 64C0C4 — 16C| |. From these expressions, we 
can derive a number of useful formulae. For the length 
of the contour line of f(x,y) at height / — h, we find 



L(h)= / |V/| p^,V/)d 2 (V/) = v^p(M (15) 



The fraction of the total sample area which lies within 
that contour is 



W(h) = / p(f) df = 



1 



(16) 



and the total Laplace curvature of / within that area is 



h +00 



K(h) 



J J A/ „(/, A/) d(A/) df = AC 2P (h) (17) 



—00 —00 



In order to fulfill the boundary condition, eq. (J™, the 
vertical position of the three-phase contact line, which 
may be symbolically written as f(dW), will vary along 
dW about an average value, < f{dW) >■ The three- 
phase contact line will thus approximately follow the con- 
tour line at /(x) =< f(dW) >, with excursions towards 
both the outside and the inside of W. These will in cases 
represent detours, sometimes shortcuts with respect to 
dW. As a reasonable approximation, we may thus use 
L(< f(dW) >) for the length of the three-phase contact 
line. Similarly, we set 



/ dx re W{h) 
Jw 



(18) 



/ A/ dx ps K(h) 
Jw 



(19) 



with h :=< f(dW) >. Inserting these expressions in 
eq. pi), we obtain 



tan 



K{h) - 2HW(h) 



(20) 



from which h can be determined. Inserting (151 and ( 17 ) 
in eq. (21 ), we arrive at 



tan# 



16C 2 _ 2HW(h) 
7T L(h) 



(21) 



B. The phase diagram 

If the adsorbed material is at liquid- vapor coexistence, 
the mean curvature of the free liquid surface, H, vanishes 
everywhere on W . In this case, eq. (21 1 is fulfilled only 
for a certain contact angle, 



8 C = arctan( a/16C 2 /7t) 



(22) 



Note that 8 C is independent of h. This at first glance 
puzzling result has its origin in a peculiar property of 
Gaussian roughness, namely the statistical independence 
of V/ and / [5H US]- I n other words, the probability 
of finding a certain slope at a given level, h, is indepen- 
dent of h. A von Neumann boundary condition such as 
eq. Q can thus be fulfilled equally well at all levels of 
Gaussian roughness. It is therefore not surprising that 
no particular value of h is singled out here. 

It is instructive to compare 9 C with 9w- For Wenzel's 
parameter r, we have 



I Vl + (V/) 2 d 2 x ps 1 + 2C 2 
Js 



0(C\) (23) 



and therefore cos 2 9w ~ 1 — 4C 2 . For C we obtain, from 
eq. Q, cos 2 6> c ps 1 — I6C2/71-. Since 16/tt > 4, we 
sec that 9 C > 6w In other words, if the liquid does 
not wet the substrate well enough to fulfill the Wenzel 
condition, is may nevertheless well intrude the rough- 
ness topography and thus form a wetting layer. This 
is indicated in Fig. [l] which shows the phase diagram 
of wetting on a surface with Gaussian roughness. States 
corresponding to liquid/vapour coexistence lie on the ver- 
tical axis. Along the bold solid line, which ends at 9w, 
the liquid surface can 'detach' from the rough substrate, 
such that an infinitely thick liquid film may form. For 
9\y < 9 < 9 C , the liquid/vapour interface needs the sup- 
port of the spikes of the roughness, to which it is attached 
by virtue of the boundary condition, eq. ([6]). 

Let us now discuss what happens as we move off coex- 
istence. We first define the parameter 



A = ^16C 2 /tt - tan 9 



(24) 
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which only depends upon the substrate topography 
(through Ca) and 9. Inserting this into eq. (21), we ob- 
tain 



2H W(h) = A L(h) 
as an alternative form of ( plj ). 

tan6> 

dry 



(25) 



tan# c 



non- 
spreading 



spreading 



H 



FIG. 1: The phase diagram of wetting and spreading on a 
surface with Gaussian roughness. 

For Gaussian roughness, it is well known that percola- 
tion of the set "P(/i) = {(x,y) \ f(x,y) < h} takes place 
at h = [2"6] . From this and eq. (15), (16), and (25) we 



can immediately derive that percolation will take place 
when 



H = A 



9l_ 
2C 



(26) 



This is indicated by the dashe d straight line in Fig. [T] 
which for 9 = ends at H c = ^/%C\I-kC$. This line cuts 
through the whole range of contact angles below 9 C . At 
higher angles, eq. ^ cannot be fulfilled and the substrate 
remains dry everywhere. 



C. Adsorption isotherms 

We can now calculate the adsorption isotherms of the 
system, i.e., the amount of liquid adsorbed at pressures 
below saturation. This is important to discuss, as almost 
no practical situation corresponds exactly to liquid/gas 
coexistence. Consider, for instance, the substrate to be 
located at a height Z above a liquid reservoir with which 
it can exchange material. H is then given by the balance 
with the hydrostatic pressure and reads 



H = 



pgz 

27 



(27) 



For the distribution of water within a soil or granular 
pile at height Z above the water table, we find that H 
grows to about 1/100 nm as Z increases to 10 m. Hence 



the typical curvatures to expect are of the right size for 
our considerations up to several meters above the water 
table. 

A more general way to look at this situation is to 
consider the vapour pressure, which is reduced at finite 
height above the liquid reservior, as well due to gravity. 
The curvature is then given by the Kelvin equation, 



H = In 



Pi 



P 2jv„ 



(28) 



where p is the partial pressure of the adsorbed liquid 
in the sorrounding gas phase, p s is its saturated vapor 
pressure, v m its molecular volume, and ks is Boltzmann's 
constant. 

From eq. (28), with the abbreviation a = kBT/jv m , 



we can express the adsorption isotherms in terms of h as 

XL(h) 



Ps 



exp 



\>W(h) 



(29) 



We would, however, like to know not the position of the 
liquid surface, h, but the total volume of adsorbed liquid. 
The latter can be easily expressed as 



V 



(h-f)p(f)df 



1 + erf 



The volume at percolation, i.e. at ft, = 



(30) 

0, is V p = 

yJCo/2-K. Combining eqs. (29) and (30), we can directly 
plot the adsorption isotherms, which are displayed in 
Fig. [2] for three different values of A/a. If A = 0, V 
remains zero for all p < p si and jumps to infinity at 

P = Ps- 
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FIG. 2: Adsorption isotherms on Gaussian roughness. The 
parameter in the family of curves is A/a. 

It is important here to appreciate that the infinite 
adsorption one obtains at coexistence has two different 
meanings for contact angles above or below 9\y. While 
for 9 < 9\y , the liquid can detach completely from the 
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substrate forming a bulk liquid phase, the liquid surface 
remains in contact with the rough substrate for 9 > 9w- 
The fact that even then the adsorption isotherms tend to 
infinity at coexistence owes to the infinite support of the 
Gaussian distribution, eq. ([9|. We will see below that 
this is a peculiarity of Gaussian roughness, and not a 
generic feature of practically encountered roughness to- 
pographies. 



IV. NON-GAUSSIAN ROUGHNESS 

As we have seen so far, it is worthwhile to study 
non-Gaussian roughness models as well. In fact, it has 
been shown that many real surfaces are distinctly non- 
Gaussian 27-30], such that the freedom in adjusting the 
correlation function is not sufficient to describe a rel- 
evantly large class of surfaces. It seems to be widely 
believed that the correlation function together with the 
height distribution of the topography are sufficient to 
characterize all physically relevant properties of a sur- 
face. Most authors even seem to believe that only the 
first four moments of the height distribution are relevant 
(including skewness and kurtosis) [2TH2U1 |3"T] . In what 
follows, we will introduce a simple roughness model which 
is general enough to describe roughness profiles with any 
lateral correlation function and height distribution, but 
is still accessible to the analysis given above. As a con- 
sequence, we will be able to derive, by purely analytic 
methods, quite general predictions about wetting, ad- 
sorption, and liquid percolation on a rough surface, which 
can be quantitatively applied to experimental data. 

Let x( x ) be a Gaussian random function, much like / 
as discussed above, but with dimensionless codomain and 
unity root mean square. Hence its height distribution is 



sured, we can derive T(/) by means of the simple formula 



q(x) = (2tt) » exp 



1 



and the correlation function, 

oo 

< x(x)x(x + r) >= 1 - D » 

in— 1 

We then set 

/(x) = S(x(x)) 



(31) 



(32) 



(33) 



where S has the dimension of a length and is a monotone, 
two times differentiable function. In this case the inverse 
of S(x), T(f) := S~ 1 (f), exists, and we have 



P(f) = T'{f) q(T(f)) 



(34) 



where the prime indicates the derivative with respect to 
the argument. 

S(x) can be directly determined from experimental to- 
pography data. If the distribution p(f) has been mea- 



T(f) = erfo 1 



2 / p(z) dz - 1 



(35) 



Note that this allows to represent any height distribution 
function p(f). The correlation function of %, and thereby 
the set of coefficients D m , is obtained from the data as < 
T(/(x))T(/(x + r)) >. Fig. [3]shows a sketch of a typical 
S(x)- While the support of S is the whole x the 
codomain is bound, because neither will there be any 
material outside the original (unworn) surface, nor will 
there be infinitely deep troughs. 



FIG. 3: Sketch of a typical S(x), as would be obtained on 
a worn surface. While the support of S is the whole \ axis, 
the codomain is bound because neither will there be any ma- 
terial outside the original (unworn) surface, nor will there be 
infinitely deep troughs, hinf indicates the inflection point of 
S. 



Since in eq. ( 33 ) we have done nothing but distorting 
the assignment of vertical positions to the plane, contour 
lines and their enclosed areas will change in level accord- 
ing to S(x)> but their topolgical properties, including 
percolation, will remain unchanged. We can therefore 
directly write down the contour length with help of the 
new quantities, 



L s (h) =^KD 2 q[T{h)] 
For the joint probability of / and V/, we obtain 



T' 2 (Vf) 2 
4C 2 



(36) 



(37) 



where the prime now denotes the derivative with respect 
to the argument. The Laplace curvature is given by A/ = 
5"'(Vx) 2 + S'Ax, which leads to intimidatingly clumsy 
expressions when inserted into multivariate analysis. We 
therefore consider here the important case when S" is 
small, such that only the second term in A/ contributes. 



G 



This is the case if 



S> 



Do 



(38) 



Wc then have 



p(/,A/) = 



exp 



< (V/) 2 > 

64£> 4 T 2 - 8D 2 TT'Af + (T'A/) 2 



2R 



(39) 

with R —\ — I6D2 |- Now we are in shape to 

express the Laplace curvature inside the wetted area. In 
complete analogy to the derivation above, we find 



K s (h) 



4Do 



-q[T(h)] 



T'(h) 

At coexistence, we have again 

L s (h)taxi6 = K s {h) 



(40) 



(41) 



In analogy to the above discussion, we define the param- 
eter 



A := 




tan 9 



(42) 



which this time does depend on h, as is sketched in Fig. [4] 
The maximum of the curve lies at hi n f, which is the 
inflection point of S(x) ( c f- Fig. [3]). 



A. The phase diagram 

The film thickness at coexistence can be derived from 
the zeros of A, of which there are either two or none. In 
the latter case, the contact angle (and thereby tan#) is 
too large for forming a liquid surface between the spikes 
and troughs which complies with the boundary condition, 
eq. ([6]). If, however, A intersects the /i-axis, the slope of 
the zeros decides upon their stability. This can be seen 
by appreciating that A may be interpreted as a deviation 
from the force balance, eq. (41 1, as required by eq. ^ . 
For the left zero, which is marked by an open circle in 



the figure, a displacement of the three-phase contact line 
would give rise to an imbalance of wetting forces driving 
it further away from the zero. The opposite is true for 
the right zero, marked by the closed circle. The latter 
is therefore stable and thus corresponds to the adsorbed 
film thickness which will develop. The film will be per- 
colated if this zero lies to the right of ho — 5(0), which 
corresponds to the mid-plane of x(x) (cf. Fig. [3]). 

If we now again consider the system off coexistence, we 
have 



2H W s (h) = A(h) L s (h) 
as the condition for h, where 



W s (h) 



(43) 



(44) 





■Mh) 


tan 9 1 


K f h \ 


h 













FIG. 4: Sketch of A(h) corresponding to S(x) as in Fig. [3] 
The film thickness at coexistence corresponds to the zero with 
negative slope (closed circle). 



A graphical solution of eq. (431 is sketched in Fig. 5] 
Again, the closed circle indicates the stable solution. The 
liquid film will be percolated if this point lies to the right 
of the dashed line at ho, but form isolated patches oth- 
erwise. From eq. (43), we see that percolation occurs 
if 



H = 



'Mho) 



(45) 



which represents again a straight line in the phase dia- 
gram as depicted in Fig. [6j For Gaussian reoughness, we 
have S( X ) = VCoX, An = C m /C , and V = I/VCq. It 
is readily checked that this leads again to eq. ( 26 ) instead 
of (??), and p4[ instead of @. 




FIG. 5: Graphic construction for solving eq. (43 1. The dashed 
line represents the l.h.s. of eq. (43 1. 



As in the case of Gaussian roughness, 6w generically 
lies below C . This can be seen from calculating 



(cos 8w) = r- 1 w 4£> 2 / q( X )S u {xW (46) 
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which follows from (37). On the other hand 



(cos 9 C 



^S-(O) 



(47) 



Since, again, 16 /ir > 4, it is clear that whenever 
argmax(S") = T(hi n f) lies close to (which it typically 
will), we have 9 C > 9\y as for the Gaussian case (cf. 
Fig.§. 



.tan 6 




FIG. 6: The wetting phase diagram for non-Gaussian rough- 
ness. The most prominent qualitative difference with respect 
to Gaussian roughness is the existence of a 'Wenzel prewetting 
line', at which the adsorbed film thickness jumps discontin- 
uously from zero (or, in fact, a molecularly thin adsorption 
layer) to a finite value. 

Inspection of Fig. [5] shows that the two points of inter- 
section, which are marked by the closed and open circles, 
will merge when the dashed and solid curves touch each 
other in a single point. This occurs at a certain curvature 
Hp (8) of the liquid surface. For H > H p , there is no liq- 
uid adsorbed, and the substrate is dry. As H p is reached, 
the average position of the liquid surface, h, jumps dis- 
continuously to the value given by the point of contact 
of the two curves. As H is further reduced, h increases 
until at coexistence it reaches a value corresponding to 
the right zero of AL$. Because of the phenomenologi- 
cal similarity of the jump in adsorbed film thickness to 
the prewetting transition encountered in standard wet- 
ting scenarios on flat substrates |16j . we hereby propose 
to term this transition 'Wenzel prewetting'. When the 
microscopic contact angle is varied, a 'Wenzel prewet- 
ting line' results, which is shown in Fig. [6] as the solid 
curve. As in the usual prewetting scenario, this line ends 
in a critical end point, when the solid and dashed curves 
in Fig. [5] intersect in only a single point. It is readily ap- 
preciated from the construction sketched in Figs. [4] and [5] 



however, that this can occur only for 9 < 0, and thus out- 
side the physically accessible parameter range. In prin- 
ciple, the Wenzel prewetting line may intersect the per- 
colation line. The latter then follows the prewetting line 
down to 9 = 0. 




Let us discuss how the position of the liquid surface 
varies along liquid/ vapour coexistence as 9 is gradually 
decreased from above 9 p . This can be directly read off 
Fig. |4j by inverting A(h) for h > hi n f, and is sketched 
in Fig. [7] The liquid film first appears through a discon- 
tinuous jump as 9 crosses the Wenzel prewetting line. 
As 9 < 9w, the liquid surface configuration which is 
bound to the surface topography through eq. ^ becomes 
metastablee (dashed curve), and the global minimum of 
the total free energy corresponds to the 'detached' liquid 
surface, or bulk liquid adsorption (bold line in Fig. [7]) . 

It is tempting to try to calculate the macroscopic con- 
tact angle, 9 macro , along the coexistence line for 9w < 
9 < 9 p . In that range, the fraction Ws(h) of the sam- 
ple is covered with liquid, while the remaining fraction, 
I — Ws(h), still exposes the uncovered rough substrate. 
The liquid surface energy of the areas covered with liquid 
is Ws(h)G(h), where G(h) =< y/l + (V.9) 2 > w is the 
total liquid surface area over W. With the help of (46) 
we readily obtain 



cos 6 n 



G{h)- / q( X ) [G(h)-4D 2 cos 9 S' 2 ( X )]d X 



T(h) 



Unfortunately, there is no straightforward way to calcu- 
late G(h). We thus content ourselves here with an upper 
bound for 9 macro , which is obtained by setting G = 1 in 
the above expression. Qualitatively, we can nevertheless 
conclude that since the jump at the prewetting line will 
directly enter in the lower boundary of the integral, it is 
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clear that this jump will as well be visible in 9 maC ro(9)- 
This is in contrast to, e.g., first order wetting, where a 
jump in film thickness is accompanied by a continuous 
variation in the contact angle [TB]. We mention again 
that 9 macro may be subject to significant contact angle 
hysteresis [2] unless long equilibration times are taken 
into account. 



B. Adsorption isotherms 

It is finally instructive to discuss the qualitative shape 
of the adsorption isotherms in this scenario, which are 
sketched in Fig. [8] The curves, which are meant to cor- 
respond to different values of £?, follow what one would 
expect for the characteristic shown in Fig. [3] The jump 
from zero film thickness to a finite value is generic and oc- 
curs for all contact angles. As 9 is increased, the height of 
the jump increases slightly, following the curvature of the 
maximum of A(/i). At the same time, the maximum film 
thickness (reached at coexistence) decreases, until it fi- 
nally comes below the percolation threshold when 9 > 9 C . 
When 9 > 9 p , A(h) has no zero anymore, and the sub- 
strate remains dry up to coexistence. 



v_ 

Vp 



1 1 1 1 — 

percolation 


A 


i ! — ■ 1 — ■ 




i 



0,0 0,2 



0,4 0,6 

Pfrs 



0,8 1,0 



FIG. 8: Sketch of a family of adsorption isotherms on non- 
Gaussian roughness. The jump occurs at what we shall term 
the 'Wenzel prewetting line'. The curves correspond, from 
bottom to top, to 9 C < < P , 9 W < 9 < 9 C , < 9 < 0w, 
and 9 = 0. The latter curve ends at finite adsorbed volume, 
but with infinite slope. 

A few more words concerning the shape of the function 
A(/i) are in order. If /(x) were Gaussian, S(x) would be 
just of the form \fCoX- I n that case, A(h) would in Fig. [4] 
be represented by a horizontal line above the abscissa. 
The solid curve in Fig. [5] would then be a Gaussian, and 
the adsorption isotherms would of course be the same as 
in Fig. [2] However, the fact that any real roughness is 
bounded, as there are neither infinitely high spikes nor 
infinitely deep troughs, entails the boundedness of the 



codomain of S(x), in contrast to the infinite codomain of 
X- As an immediate consequence, the derivative of T{h) 
must finally diverge at the boundaries of its (finite!) sup- 
port, which necessarily leads to A bending down onto the 
dashed line below the abscissa in Fig. [4j This leads not 
only naturally to a finite h at coexistence (H — 0, cf. 
Fig. [5]), but also to the Wenzel prewetting jump in h(H) 
farther away from coexistence, when the dashed curve in 
Fig. [5] just touches the solid curve. This reveals that the 
shape of the adsorption isotherms we derived for Gaus- 
sian roughness above is qualitatively different from what 
should be expected for real roughness. In fact, it misses 
the whole prewetting scenario, which turned out above 
to be a generic feature. Once again, Gaussian roughness 
reveals itself as a special case, which is mathematically 
convenient but may be misleading when it comes to mak- 
ing quantitative predictions. 



V. CONCLUSIONS 

In conclusion, an analytic theory was presented which 
allows to calculate the wetting phase diagram, adsorption 
isotherms, and percolation threshold of the adsorbed liq- 
uid film for isotropic, randomly rough substrates with 
arbitrary lateral correlation function and height distri- 
bution. The results are found to depend only upon a 
few key parameters, which can be clearly identified and 
derived from experimental sample profile data. We have 
seen that wetting 'physical' roughness displays a num- 
ber of features which are not present for exactly Gaus- 
sian roughness, such as a prewetting transition occur- 
ring well before the Wenzel angle is reached. This could 
be traced down to subtle properties of Gaussian random 
functions, which reveal their unphysical nature only at 
second glance. 

Since for most quantities of interest we could come 
up with closed analytic expressions, these results may 
be particularly useful for practical applications. The 
range of validity of the present theory extends from a 
few nanomeres up to roughly a millimeter, well below 
the capillary length of the liquid. The field of such appli- 
cations is vast, including almost all situations in which a 
liquid comes into contact with a naturally rough surface. 
In particular, ramifications of wetting phase transitions, 
which inherently involve small contact angles, are to be 
expected and can now be accounted for in closed form. 

Given the potential relevance of the results presented 
here, it will be worthwhile to work on relaxing the five 
major approximations we have used: 

1. We have assumed the substrate to be chemically 
homogeneous. 

2. We have assumed shallow profiles; eqs. Q and d5]). 

3. We have neglected higher correlations in the shape 
of dW; eqs. (18), (19), and the paragraph before. 
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4. We have assumed the curvature of the roughness 
characteristic S(x) to be small; eq. (38). 



5. We have assumed isotropy of the roughness; 
eq. flit. 



The last two are probably the simplest to tackle, while 
the first one appears as the most difficult to overcome. It 
should finally be noted that in experiments, one has to 
reckon with quite long relaxation times for the measured 
quantities, because at all levels of the roughness there 



are saddle points [STJ H3] , which act as effective pinning 
sites for h. Equilibration will nevertheless proceed within 
manageable time, either via the vapour phase or via the 
molecularly thin adsorbed film [S]. 

Inspiring discussions with Daniel Tartakovsky, 
Siegfried Dietrich, Martin Brinkmann, Jiirgen Vollmcr, 
Sabine Klapp, and Daniela Fliegner are gratefully 
acknowledged. The author furthermore acknowledges 
generous support form BP International. 
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